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I. INTRODUCTION 

The elastic KK scattering is one of the simplest re- 
actions with strange quark, and it allows for an explicit 
exploration of the three-flavor hadronic structure. The 
measurement of KK scattering is very useful for our 
study of the chiral symmetry breaking of quantum chro- 
modynamics (QCD). In the isospin limit, the KK system 
have two isospin eigenstates, namely, 1=1 and 

During the last decade, the experimental measure- 
ments of the KK scattering have been carried out at 
several cooler-synchrotron COSY facilities Because 
of the high momentum resolution of COSY, the thresh- 
olds from neutral and charged kaons are well separated 
and can be studied independently, and the values of to- 
tal and differential cross sections are now available for 
a variety of reactions [U-Q, where the prompt K + K~ 
spectrum is described by a four-body phase-space distri- 
bution, modified by the various final state interaction [2- 
§]. These results show a clear evidence for the smooth 
K + K~ background. 

It should be stressed that, based on the low en- 
ergy K + K~ invariant mass distributions and gener- 
alized Dalitz plot analysis, Silarski ct al. first esti- 
mated the scattering length for the K + K~ interaction 
to be \Re(a K +K-)\ = 0.5^0.5 fm and Im{a K + K -) = 
3.0 ± 3.0 fm Q. Later, they improve these results as 
]Re{a K+K -)\ = 0.2;°;* fm, Im{a K+K -) = 0.4±g;2 fm 
[0| . Since K + K~ system is a mixture of the isospin 1 = 
and 1, in those analyses we can not distinguish between 
the 1 = and 1 channels. 

The determination of KK scattering from QCD is very 
difficult since it is essentially a non-perturbative prob- 
lem. However, some theoretical efforts are still taken to 
study the KK scattering [lol - [T2l |. But, if the scattering 
hadrons contain strange quarks, the chiral perturbation 
theory (%PT) predictions usually suffer from considerable 
corrections because of the resonances / and a , and one 
has to apply chiral unitary theory, sometimes called uni- 
tarized chiral perturbation theory (i~3l - [l6j . 

The most feasible way to extract KK scattering length 
nonperturbativcly from first principles is using lattice 



QCD. And it offers an another important consistent 
check of the validity of xPT with the inclusion of the 
strange quarks. Although there are some exploratory 
lattice QCD investigations of the meson anti-meson po- 
tential ( including KK potential ) in Refs. [13, El: until 
now, no lattice QCD study about its scattering length 
has been reported, mainly because it is extremely dif- 
ficult to reliably calculate the rectangular and discon- 
nected diagrams. Motivated by the first lattice study of 
the s-wave K + K + scattering length in the 1=1 channel, 
which was explored by NPLQCD Collaboration in fully- 
dynamical lattice QCD with domain- wall valence [l9| and 
the value of mK+ a K+K+ was found to be —0.352(16), we 
will explore the KK scattering length from lattice QCD. 
This is also encoura ged by our reliable extractions of the 
irK [2(|[2l| and irn [22j| scattering lengths. 

Luscher [23|, [24[ established the basic formulae for the 
calculation of the scattering length using lattice QCD, 
which is valid for the elastic scattering below inelas- 
tic thresholds. Below kinematic thresholds, the scat- 
tering length of two hadrons is connected to the energy 
phase shift of two-hadron state enclosed in a torus. This 
method paved the way for our lattice calculations of KK 
scattering. We should bear in mind that, above the in- 
clastic threshold, a tower of resonances emerges, which 
suggest the opening of other channels. 

It is well-known that the appearance of the fo(I = 0) 
and ao(I = 1), is due to the introduction of the KK scat- 
tering channel [l|. They decay mainly into tttt and TVn, 
respectively. Both have the masses around 980 MeV [25| ■ 
Since the central mass values fall almost exactly at the 
KK threshold, the strong coupling to this channel dis- 
torts significantly the upper parts of the mass spectra. 

As presented later, for isospin 1 = channel, we need 
calculating vacuum diagram, which is extremely difficult 
to reliably measure. Thus, we here will preliminarily re- 
port our lattice simulations on the KK scattering length 
in the 1=1 channel. In the presence of the strange 
quark exchange channel, the KK scattering length have 
imaginary contribution. For instance, the experimental 
determinations of the scattering lengths are complex as 
dictated by the various decay channel thresholds 043- 
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The S matrix has the structure 



S = 



j(l _ ^2)1/2 e i(5 K ^+<5^) 



where 8 K j^ and are the phase shifts for the elastic 
KK — > ATAT and nrj — > 7r?/ processes in isospin 7 = 1 
channel and rj is the inelasticity. We note that, using 
(S)u and (S) 22 we can determine 77, t^g* and 5^. 

We will use the MILC gauge configurations gener- 
ated with the Nf =2 + 1 flavors of the Asqtad im- 
proved [H,[23] staggered dynamical sea quarks [28|,|2£| to 
calculate the elastic KK phase shift 8 K -^ and then eval- 
uate the s-wave KK scattering length for isospin 1=1 
channel in the absence of resonant states. We adopt the 
technique introduced in Refs. [30l. [3l|. namely, the mov- 
ing wall sources without gauge fixing, to obtain the re- 
liable accuracy, and observe a clear signal of attraction, 
which is consistent with pioneering lattice studies on KK 
potential in Refs. [13, Ea|- Moreover, we extrapolate the 
KK scattering length to the physical point using the con- 
tinuum three-flavor xPT form at the next-to-leading or- 
der (NLO), which are provided in [B] with the help of 
Zhi-Hui Guo, and it can be directly built from the its 
scattering amplitudes in Ref. [HI], 

This article is organized as follows. In Sect. [Til wc de- 
scribe the formalism for the calculation of the KK scat- 
tering lengths using the Liischer's formula [23I . |2~H ] and 
our computational technique for KK four-point func- 
tions. We present our lattice results in Sect. IIII1 and 
arrive at our conclusion and outlook in Sect. IIVI 



II. METHOD 

Here we follow the original notations and conventions 
in Refs. [30l - l34j to review the required formulae for cal- 
culating KK scattering lengths. In the Asqtad-improved 
staggered dynamical fermion formalism, let us study the 
KK scattering of one Goldstone kaon and one Goldstone 
anti-kaon. Using operators Ok{xi), Ok{%2) for kaons at 
lattice points X\, 22, and operators 0^(23), 0^(24) for 
anti-kaons at lattice points 23, 24, respectively, we can 
express the KK four-point correlation function as 

C K k(x4,x 3 ,x 2 ,x 1 ) = (O w (x4,)O w (x s )0] ( (x2)Ol c (x 1 )), 

(1) 

where the kaon and anti-kaon interpolating field opera- 
tors are denoted by 



o K+ {t) 

<D K o(t) 

O^o (t) 
K -(t) 



^s(x,t)7 5 w(x,t), 

X 

^s(x,i)7 5 d(x,t), 

X 

-^2 d(x,f)7 5 s(x,i), 

X 

^u(x,t)7 5 .s(x,t). 



(2) 



After summing over the spatial coordinates, we achieve 
the KK four-point correlation function, 

C K K(t4,t 3 ,t2,ti) = ^2^2^2^2C K k(x4,X 3 ,X 2 ,X 1 ), 

Xi X 2 X 3 X 4 

(3) 

where x\ = (xi,fi) (likewise for 22, 23, and 24), and t 
represents the time difference, namely, t = t 3 — t\. We 
build KK operators for two isospin eigenstates as [l| 1 

\KK(I = 0)) = Jl{\K+K-) \K°K^)}, 

\KK{I=l)) = ^-{\K+K-)+ \K°W)). 

Then, the scattering amplitudes M for KK scattering in 
the 1 = 1 and states are given by 



M(I=1) = -(K+K-\S\K + K- 
+ ^(K°K 1 >\S\K + K-] 



h^(K + K~\S\K°K°) 
hK°K°\S\K°K°),(4:) 



M(I=0) = ^(K+K-\S\K+K-) - ^(K + K-\S\K°K°) 

-^(K°K°\S\K+K-) + ^(K°K°\S\K KV).(5) 

We substitute Eq. ((2J into Eq. ([I]) to obtain the quark 
diagrams for the KK scattering in the 1=1 channel. 
The concrete calculations are given in [X] for reference. 
At last, we obtain 12 different diagrams for the 1=1 
channel, which correspond to Eqs. (|A1|) - (|A12[) . and are 
also shown in Figs. [TJ and [21 
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FIG. 1: The numbers below the quark line diagrams cor- 
respond to the terms (|A1|) - (|A6|) given in Appendix A. (A) 
and (B) correspond to the first and second terms in Eq. Q, 
respectively. 



1 Our phase conventions for pscudoscalar mesons arc different 
from those in Ref. Q. 
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FIG. 2: The numbers below the quark line diagrams corre- 
spond to the terms (|A7|) - (|A12|) given in Appendix A. (C) 
and (D) correspond to the third and fourth terms in Eq. 
respectively. 

In the isospin limit, we can classify the 12 quark line di- 
agrams into four independent groups. Diagram 1 in Fig.Q] 
and diagram 9 in Fig. [5] arc categorized into Group 1. 
Likewise, diagrams 3 and 10 are categorized into Group 
2; Nos. 2,5,7, and 11 in Group 3; Nos. 4,6,8, and 12 
in Group 4; respectively. By inspecting Eq. (j4|) and 
Eqs. (|A1[) - (|A12[) . the coefficients of these groups are 
given as 



Group 1: 
Group 2: 
Group 3: 

Group 4: 



1 1 

2 + 2 =1:D < 
§(-l) + |(-D 



— —1 : R u , 



(6) 



(-1) + 
(-1) + 



where, R u , R s are two types of the rectangular quark dia- 
grams and the subscript u and s distinguish two of them, 
since for R u rectangular quark diagram, there are two u 
quark lines on the time direction, likewise for R s . At 
last, only two diagrams still contribute to KK scattering 
amplitudes for the isospin 1=1 channel. For the 1 = 
case, we can easily perform the same procedures. Finally, 
both channels can be written in terms of the diagrams 
D, R^, R s and V, namely [13,1111 

M(I = 1) = D — R u , 

M(I = 0) = D - Ru - 2R S + 2V. (7) 

To avert the Fierz rearrangement of quark lines (30l . 
Hl|, we select ti = 0, t 2 = 1, t 3 = t, and f 4 = t + 1. 



Then, we can build KK operators for isospin 1=1 and 
eigenstates as [l| 



^(£) = -±= {K+(t)K-(t + l)-K°(t)K°(t + 1)} , 
Kit®=\ {K + (t)K-(t + l)+K°(t)K°(t + 1)} ,(8) 
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FIG. 3: Diagrams contributing to KK four-point functions. 
Short bars stand for wall sources. Open circles are sinks for 
local kaon or anti-kaon operators. The blue and red lines 
represent the u/d and strange quark lines, respectively. 

In the isospin limit, four quark line diagrams con- 
tribute to the KK scattering amplitudes; we plot them 
in Fig. [31 and label them as direct (D), rectangular 
(R u ), rectangular (R s ), and vacuum (V) diagrams, re- 
spectively It is pretty easy to evaluate the direct di- 
agram, while the reliable evaluation of the rectangular 
(R u ,Rs) and vacuum (V) diagrams is extremely diffi- 
cult [13, HU . We settle it through the moving wall source 
technique introduced in Rcfs. [30|,[3l[, namely, each prop- 
agator, which corresponds to a moving wall source at 
t = (),■■ ■ , T — 1, is denoted by 



',(x,t)> 



< t < T- 1, 



where D defines the quark matrix for the quark action. 
D, Ru, R a and V are schematically displayed in Fig. 02 
and we can also write them in terms of the quark prop- 
agators G, namely, 

C D (U,ts,t2,h) = ^^<Tr[G^ )t (x3,t 3 )G^(x3,i 3 )] 



Tr[Gi: )t (x 4 ,t 4 )G^(x 4 ,i 4 ) 



0), 



C Ru (t 4 ,t 3 ,t 2 ,h) = ^^(Tr[G|; s)t (x 2 ,t 2 )G| 4 u) (x 2 ,t 2 ) 

Cr s (*4j *3j *2) *l) 



X2 X3 

xGi; )t (x 3 ,t3)G(; ) (x3,t 3 )]>, 

^^(Tr[G^ )t (x 2 ,f 2 )GW(x 2 ,t 2 ) 



xGi 4 U)t (x3,t 3 )G^(x3,t3)]>, 

Y,^{mGif(x 2 ,t 2 )Gi:\x 2 ,h)) 

xTr[G| 4 u)t (x3,t3)G(: ) (x3,i 3 )]>- 
(Tr[Gif(x 2 ,t 2 )G£\x 2 ,t 2 )}) 

x(Tr[Gi 4 u)t (x 3 ,t 3 )G(: ) (x3,i3)])},(9) 



Gy(i 4 , t$, t 2 , ti) 
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where we utilize the hermiticity properties of G to re- 
move the 7 5 factors. The vacuum diagram here include 
a vacuum subtraction term. 

According to the discussions in Refs. [1q|, HH], the KK 
rectangular and vacuum diagrams in Fig. [3] produce the 
gauge-variant noise, and we will perform the gauge field 
average without gauge fixing to nicely suppress this noise. 
We can construct physical KK correlators with certain 
isospin by integrating four types of the propagators. In 
the isospin limit, we can express the KK correlator for 
isospin 1=1 and channels by 



KK^ 



KK 
C I=1 



<0*#(t)|0g(o)) 



D 
D 



Ru, 



2iL 



2V, 
(10) 



where operator O k -jt creates a KK state with total 
isospin /. Since, in this paper, we preliminarily report 
our lattice simulations on KK scattering length in the 
1 = 1 channel, in the following discussions, we will re- 
move the superscript / for the corresponding quantities. 

There are some different methods to parameterize the 
low-momentum behavior of the scattering amplitude. To 
calculate the KK scattering lengths on the lattice, we 
adopt the Liischer formula, and it is straightforward to 
use the standard effective range expansion for the KK 
scattering phase shift, namely, 

k cot S (k) = - + lrk 2 + 0(k 4 ), (11) 
a 2 

where So(k) is s-wave scattering phase shift, a is the scat- 
tering length, r is the effective range, and k is the magni- 
tude of the center-of-mass scattering momentum related 
to the total energy o f the KK system in a cubic box of 
size L by E K -^ = 2yJm z K + fc 2 , where vtlk is kaon mass. 
This expansion is only valid at low energy, unfortunately, 
takes inelasticities not into consideration. 

If we approximate the KK element of the S matrix by 
its dominant pole contribution, the s-wave KK scatter- 
ing length in the continuum is denoted by employing the 
standard Liischer formula 



ao 



lim 

fe^O 



tan^o(fc) 



where the scattering is purely elastic below inelastic 
thresholds. We should bear in mind that the trunca- 
tion of the effective range r in Eq. (|llj) serves as a source 
of the systematic uncertainty, which appears as 0(1/ L 6 ) 
and in this paper we ignore it. 6o(k) is the s-wave scatter- 
ing phase shift, which can be evaluated by the Liischer's 
finite size formula K23l \2$ , 



tani5o(fc) 



4tt 



■z. 



mi 



1, 



k 2 



k J ttL \'' (2tt/L) 2 i 

where the zeta function Zqq(1] q 2 ) is denoted by 



(12) 



(13) 



here q = kLj (2tt) , and in this work the zeta function 
Z 00 (l; q 2 ) is calculated by the method in Ref. [35| . 

The energy E K ^ of the KK system can be obtained 
from the KK four-point correlator. At large t these cor- 
relators behave as [39, K37f 



C 



KK 



E Klz \t 



(-l) 4 ^cosh 



E 'kkV-\ 



•(14) 



where E K j^ is the energy of the lightest KK state. The 
second term alternating in sign is a peculiarity of the 
staggered scheme [H, H?} . 

In our concrete calculation we also measure the energy 



shift 6E = E 



R x {t) 



— 2itlk from the ratio 
(0,1, M + l) 



KK 



c K (Q,t)c K {i,t + iy 



X=D&ndR u , (15) 



where Cr- is the kaon two-point correlator. Considering 
equation (|10[) , we can write the amplitudes which project 
out the 1 = 1 isospin cigcnstate as 



R(t) = R u (t)-R R -{t) 



(16) 



We should stress that, as compared to the contribu- 
tions of the Nambu-Goldstone kaon and anti-kaon, these 
of the non-Nambu-Goldstonc kaons and anti-kaons in the 
intermediate states are exponentially suppres sed for large 
times due to their heavier masses [30l . l3ll . |33| . In the cur- 
rent study we assume that the KK interpolator does not 
couple significantly to other KK tastes, and ignore this 
systematic error. 



III. SIMULATION RESULTS 

We calculated KK correlator on MILC lattice en- 
semble of 200 16 3 x 48 gauge configurations with bare 
gauge coupling 10/g 2 = 6.572 and bare quark masses 
am u d/am s = 0.0097/0.0484. Its physical volume is 
about 2.5 fm. The inverse lattice spacing a -1 = 1.358^3 
GeV [H, [29[ ( about 0.15 fm ). The dynamical strange 
quark mass is quite close to its physical value [H, [29| . 

The necessary matrix element of inverse fermion ma- 
trix are computed using the standard conjugate gradient 
method. We calculate the KK correlators on all the time 
slices for both source and sink. After averaging over all its 
possible values, the statistics are dramatically enhanced 
because we can place the kaon and anti-kaon sources at 
all time slices. 

Using the same lattice gauge configurations, we cal- 
culate the KK four-point correlation functions with six 
light u valence quarks, namely, am x = 0.0097, 0.01067, 
0.01261, 0.01358, 0.01455 and 0.0194, where we adopt 
MILC convention: m x is the valence u/d quark mass. 
The strange sea quark mass is chosen at its physical num- 
ber di. 
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In Fig. [4] the individual ratios, R x (X = D and R u ) 
are shown as the functions of t for am x = 0.0097. The 
values of the direct amplitude R D is quite close to unity, 
suggesting the weak interaction in this channel. After 
a beginning increase up to t ~ 4, the rectangular am- 
plitude shows an approximately linear decrease up until 
t ~ 20, indicating an attractive force between the kaon 
and anti-kaon in this channel. These features are what 
we expected from the theoretical predictions [33j | . The 
systematically oscillating characteristics in the large time 
region is also distinctly noticed, which is a speciality of 
the staggered scheme j3(| H?} . 
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FIG. 4: Individual amplitude ratios R x (t) for KK four- 
point function measured by the moving wall source without 
gauge fixing as functions of t: Direct diagram shifted by 0.8 
(diamonds) and rectangular (octagons) diagrams. 

In our previous works [2(1 [H, [39| , we have measured 
the pion and kaon point-to-point correlators, and calcu- 
lated the pion mass m w , kaon mass mx, and pion de- 
cay constants f w , which are summarized in Table I of 
Ref. [20(. In this work, we will directly quote these val- 
ues. 

According to the arguments in Ref. |4(j, the ratios 
for the rectangular R u diagram have errors, which in- 
crease exponentially as e™ 1 * 1 from large time separation 2 . 
The magnitude of the errors is in quantitatively agree- 
ment with this expectation as displayed in Fig. [SJ Fit- 
ting the errors SR Ru (t) by a single exponential fit ansatz 
SR R "(t) ~ e MRt , we can obtain the corresponding fitting 
values of hr. The fitted values of fj,R in lattice units and 



2 The errors for vacuum amplitudes should be roughly indepen- 
dent of t, and grows exponentially as e 2mK * in the ratio Ry . 
While for the rectangular diagram R s , it increases as e maSt , 
where m S s is a fictitious meson with two valence quarks with 
mass about m s . Thus, the reliable calculations of these terms 
are beyond the scope of this paper since it requires a substantial 
amount of computing resources. 
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FIG. 5: The errors of the ratios R Ru (t) as the functions of t. 
Solid lines are single exponential fits. 



its fitting ranges are summarized in Table HI From Ta- 
ble HI we can note that the fitted values of /.ir can be 
compared with the corresponding pion masses listed 
in Table I of Ref. (20|. This demonstrates, on the other 
side, that the technique of the moving wall source with- 
out gauge fixing used in the current work for the KK 
scattering is practically feasible. 



TABLE I: Summary of the fitted values of fj,R in lattice units. 
The second block shows the fitted values for Column three 
shows the time range for the chosen fit. 

am x a^in Range 



0.0097 


0.2768 


8-14 


0.01067 


0.2825 


8-14 


0.01261 


0.3047 


8-14 


0.01358 


0.3134 


8-14 


0.01455 


0.3183 


8-14 


0.0194 


0.3518 


8-14 



As we practiced in Ref. [20(, we understood that the 
correctly extracting these energies is very important to 
our ultimate results, and in our concrete simulation, 
they were chosen by searching for the combination of a 
"plateau" in the energy as the function of the minimum 
distance, and a good chi-square (namely, x 2 ). In Fig. [6] 
we plot the KK correlation function for am x = 0.0097, 
where we can compare the fitted functional form with the 
lattice data. The fitted values of the energies aE, fitting- 
range and fitting quality are listed in Table HT1 

We now can plug these energies in Table [TT| into 
Eq. pT|) to achieve the scattering lengths. The center- 
of-mass scattering momentum k 2 in GeV and the scat- 




the its scattering amplitudes in Rcf. |32j, namely, 



FIG. 6: The KK correlator calculated with the moving wall 
source without gauge fixing for am x — 0.0097. Solid line is 
the fit for 6 < t < 12 using the fitting model in Eq. (fT4")) . 



TABLE II: Summary of the lattice results for the fitted values 
of the energy. The third block shows its fit range, and the 
fourth block gives its fit quality x 2 /dof. 



am x 


aE 


Range 


X 2 /dof 


0.00970 


0.7862(6) 


6-12 


4.29/3 


0.01067 


0.7931(6) 


6-12 


4.81/3 


0.01261 


0.8071(6) 


6-12 


4.96/3 


0.01358 


0.8139(6) 


6-12 


5.01/3 


0.01455 


0.8207(6) 


6-12 


5.12/3 


0.01940 


0.8538(5) 


6-12 


5.58/3 



K 



16^/2 



h^5(M) 

J IT 



32m 



(17) 



(2 KK 

J 7T 



where wc substituted the pion mass m ff , kaon mass m 
and pion decay constant f„ listed in Table I of Ref. [2 
L 5 (n) and L^Qi) = 2L 1 + 2i 2 + J L 3 -2i 4 - |L 5 +2L 6 4 

and ex 
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Lg are low-energy constants denoted in Ref. 
plicitly dependent on the chiral scale fi. The X^^(m) is 
the known functions at NLO including chiral logarithm 
terms, see Eq. (|B5j) for details. 




m z K [deV a ] 



tcring lengths are tabulated in Table [TTTT Here we utilize 
the kaon masses given in Table I of Ref. [2(| . The errors 
of the scattering momentum k and the scattering lengths 
are estimated from the statistic errors of the energies E 
and kaon masses tuk- 



TABLE III: Summary of lattice results for 
lengths. Column two gives the center-of-mass 
mentum k 2 in GeV, and the third block shows 
times the scattering lengths. 

fc 2 [GeV 2 



the scattering 
scattering mo- 
the kaon mass 



m K a Kl y 



0.00970 


-0.00444(54) 


0.363(52) 


0.01067 


-0.00446(53) 


0.368(52) 


0.01261 


-0.00451(52) 


0.378(52) 


0.01358 


-0.00483(53) 


0.414(54) 


0.01455 


-0.00471(54) 


0.406(55) 


0.01940 


-0.00488(54) 


0.441(58) 



In Fig. the s-wave KK scattering lengths mxa 



i=i 

KK 

are displayed as a function of m K . In this work, we use 
pion masses m„ = 330 — 466 MeV, and need to extrapo- 
late the KK scattering lengths toward the physical point. 
For this end, in [B]we provide the continuum SU(3) xPT 
form at NLO for a]^, which can be directly built from 



FIG. 7: m^-dependence of the KK scattering lengths 
mKd K -ji for the I — 1 channel. The cyan diamond point 
indicate its physical values. 

To enhance the xPT fit, we include all the lattice 
simulation data of the KK scattering lengths. The fit- 
ting results of the KK scattering lengths m^a^i are 
plotted by the dotted lines as a function of m K in 
Fig. [7] The chirally extrapolated KK scattering length 
rriKa 1 ^ = 0.211(33), the cyan diamond points in Fig. [7] 
demonstrates this value. The fit parameters L. k k,Lc >1 
and the scattering lengths rriKa 1 ^^ at the physical points 
(namely, m n = 0.140 GeV, m K = 0.494 GeV) ^ are 
also summarized in Table uVl and the chiral scale /j, is cho- 
sen as the physical r\ mass, namely, \i = 0.548 GeV (2f| 
as it is done in Ref. [201 ■ 



TABLE IV: The fitted s-wave scattering lengths uikcikK 
at the physical point (m n = 0.140 GeV, m K = 0.494 GeV). 
The chiral scale fi is taken as the physical 77 mass. 



X 2 /dof 


in 3 • T I=1 


10 3 • u 




0.656/4 


-2.41 ± 1.10 


4.90 ± 4.35 


0.211 ±0.033 



From Fig. wc can observe that our lattice simula- 
tion results of the scattering lengths have a large error, 
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and are in reasonable agreement with the SU(3) x?T at 
NLO. The fitted values of the L 5 and tuk have large 
statistical errors, which reflect the big errors of our ex- 
tracted lattice data of mfca 1 ^,. Thus we can not claim 
they have physical meanings. 

In this work, we only consider the statistical errors. 
The possible sources of the systematic errors on the ex- 
trapolated value of m if mainly contains two parts. 
First, according to aforementioned discussions, when ex- 
tracting rriK we ignore two major systematic errors: 
the truncation of the effective range r and the contribu- 
tions of the non-Nambu-Goldstonc kaons and anti-kaons 
in the intermediate states. These systematic errors of the 
extracting rriKa 1 -^ should propagate through the chiral 
extrapolation. Second, the extrapolation to the physical 
point needs the experimental value for m w , ttlk and f w . 
The experimental error on these quantities brings a quite 
small systematic error as compared to the corresponding 
statistical error and therefore is neglected. 



IV. SUMMARY AND OUTLOOK 

In this work we have carried out lattice study of KK 
scattering lengths, and performed a concrete lattice cal- 
culation of KK scattering lengths for isospin / = 1 chan- 
nel, where rectangular diagram plays a vital role, for the 
MILC medium coarse (a = 0.15 fm) lattice ensemble with 
the Asqtad improved the staggered d yna mical sea quarks. 
We employed the technique in Refs. |30l.l3~fl] (namely, the 
moving wall sources without gauge fixing) to reliably cal- 
culate KK four-point correlation function, and observed 
a clear signal of the attraction for isospin 1 = 1 chan- 
nel, which is in well accordance with the pioneering lat- 
tice studies on KK potential in Refs. [13, HI]. Extrap- 
olating the lattice data of the s-wave scattering lengths 
to the physical point, we achieved the scattering length 
mxdjjj = 0.211(33) directly from lattice simulations. 

A good signal can be seen for long time separation 
range in the rectangular R u diagram of A'A' scattering. 
We can further reduce the noise by a larger lattice volume 
or smaller pion mass, and hence obtain better results for 
the scattering length. Moreover, as pointed above, if we 
want to obtain the good signals of the rectangular R s 
and vacuum diagrams, we should choose lattice ensemble 
with suitable strange quark. 

We should understand that the investigation on KK 
scattering at the 1=1 channel is just first step in the 
study of the A'A' scattering. Since the A" + A"~ system is 
a mixture of the isospin 1 = and 1, and now there are 
enough experimental results about its scattering lengths, 
the lattice study on A^A" scattering length at the 1 = 
channel is highly desired. We are planning a series of 
lattice studies on it. 
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Appendix A: I = 1 KK scattering amplitudes in 
terms of quark propagators 

(0 K+ (x 4 )O k - (x 3 )0^ K+ {x 2 )0^ (xi)) = 
Tr (g(' u >(x 2 ,x 4 )7 5 G( s >(x4,x 2 )7 5 ) x 

Tr {G^(x 3 , Xl ) l5 G^(xux 3 ) l5 ) (Al) 

- Tr(G< s) (x 2 , x 4 )7 5 GW (X4,x 3 )75 x 
G^(x 3 ,x 1 ) l5 G^(x 1 ,x 2 ) l5 ) (A2) 

- Tr (g (u) (x u x 3 ) 75 G (s) (x 3 , x 4 ) 75 x 
G (u) (x4,X2)75G (s) ( a;2 ,x 1 ) 75 ) (A3) 

+ Tr (GM (a* , x*)l$G^ (x 4 , 23)75 x ) 

Tr [g^\x 2 ,x 1 ) 1z G { - s \x 1 ,x 2 ) 1 ^ . (A4) 

(O K+ (x 4 )O K -(x 3 )O^ K0 (x 2 )O^ t (x 1 )) = 
~ Tr(G^\x 3 , Xi ) lh G^{x 4 ,x 2 ) lh ) x 

Tr^ d \x 2 ,x 1 ) l5 G^(x 1 ,x 3 ) l5 ) (A5) 

- Tr(G (s >(x3,X4)75G^(x4,X3)7 5 ) x 

Tr (G^(x 27 x 1 ) l5 G^(x 1 ,x 2 ) l5 ) . (A6) 

(0 KO (x 4 )0 1? o(x 3 )Ol (+ (x2)0^.(x 1 )) = 

- Tr(G ( ^3,Z4)75G (s) (x 4 ,x 2 ) 75 ) x 

Tr^^XihsG^^xghs) (A7) 

- Tr(G (d) (x 3 ,2 4 )75G (s) (x4, 23)75) x 

Tr (g^(x 2 ,x 1 )7 5 G( s )(x 1 ,x 2 ) 75 ) . (A8) 

(O k o (z 4 )<V 03)0^0 {x2)Oy (xi)) - 
Tr(GW(x 2 ,.T4) 75 G( s )(x4,x 2 )75) x 
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Ti(G^{x 3 ,x 1 ) l5 G^{x 1 ,x 3 ) l5j 

- Ty(g^(x 2 ,x 4 ) 15 G^(x 4 ,x 3 ) 15 x 

(Z1,X2)75 ) 

- T r (G^(x 3 ,x 4 h 5 G^(x 4 ,x 2 h 5 x 

(xi, 0:2)75 

+ Tr(G ( ^(x 3 , a; 4)75G (d) ^4,X3)75) x 
Ty(g (s) (x 2 ,x 1 ) 15 G w (an ,a: 2 )V 



(A9) 
(A10) 
(All) 
(A12) 



Appendix B: The analytic expression of s-wave i-fT-sT 
scattering length for the isospin I — 1 channel 

In the isospine limit, both isospin amplitudes in KK — > 
KK can be described by two independent amplitudes T c h 
(the amplitude for the processes K + K~ — > K + K~) and 

T neu (that for K°K° -> if+if-) [U, HH [H]. Z. H. Guo 
and J. A. Oiler derived two isospin amplitudes only in 
terms of T neu I32L Here we follow the notations and 
conventions in [33, since the the chiral scale fi depen- 
dence in amplitude is canceled by the loops and the low 
energy constants (LECs) [32| . 

These amplitudes can be decomposed into partial 
waves U(s) according to 

T(s,t,u) = 16tt^(2Z + l)tt(s)Pi(cos6), 



where I is total angular momentum, 9 denotes scattering 
angle in the center-of-mass system. In the elastic region, 
the partial wave amplitude t;(s) can be parameterized by 
real phase shifts 61 (s), 



And its real part can be expanded at threshold (namely, 
s = 4m\-, t = 0, u = 0) in terms of the scattering lengths 
(af =1 ) and effective ranges (&f =1 ), 



f( fl ) = ^^{of= 1 + 6f=V+0(3 4 )}> 



Ret 



for the center-of-mass three-momentum of the kaons q. 
The s-wave scattering length is linked to the real part of 
the amplitude at threshold by 



ReT/^ 1 



lQnrriKCL 



1=1 

KK 



0(q 2 



we first expand the scattering amplitude at threshold, 
namely 3 , 



T(s,t,u) = 



2m| 

J n 



32m|m 2 

71 L 5(P) 

Jit 



64m 



fA KK 

J IT 

.4 



(3m 2 +m 2 ) +lUm 4 K 3mi , mi 
■ — log — 



2887r 2 /^ 
1 



4327r 2 / 4 (m 2 " m D 
-ml - 228m 2 ,m 2 K 
1 



W 4 

,2 inn a 4 



100m 2 K ml 



o 2 2\ 

3m „m w ) 



43207r 2 /^(m2 - m 2 ) 
+377^(716/74- + 73m 2 ,) 

+m 2 (167m 4 - 3760tt4 
+</r-a-( s 



15m° - 441m 



log' 



n 2 
11 
,2 



a 2 \ 4:111 K 

Am 2 K )- 



+Jkt 1 {s = Am 2 K ) 



(37 



f 4 



28m? K 



54/ 4 



(Bl) 



where L^(fj,) = 2L X + 2L 2 + L 3 - 2L A -_±_L 5 + 2L 6 + L 



are low-energy constants defined in Ref. [41| at the chiral 
symmetry breaking scale fi. The loop function J pq is 
given in Ref. [4l| , namely, 



Jnri(s = Am 2 K ) 



32tt 2 



4m 2 K 



4m 



In 



(4m 



2 



(X 



v , 
,2\2 ' 



(4m|- - v) 2 - (m 2 - m 2 ) 2 



m± 



(B2) 



with ^ = 



Am 2 K 



Am 2 K 



(m„ + m v ) 

. By plugging the leading order relation m 2 = (Am 2 K 
m 2 )/3 (32( 1 we can simplify the above formula as 



T(s,t,u) = 



2m 2 K 
p 

J 7T 



Jit 

o 4 2 

7J lo S— § 



64m 



f 4 



KL K-kM 



9tT 2 / 4 87T 2 /4 

m 2 m 2 f (2m 2 f + 5m 2 ) mj. 
167r 2 / 4 (m|- - m 2 ) /x 2 

56mi 4ml 
log 



9m 4 m|- 



Wm^m 2 



144^ 2 /4( r 



A' 



+Jifjf(s = -im 2 K )^f- 



m 2 ) 



3/j 2 



The one-loop order analytic expressions of KK scat- 
terin g a mplitudes in the isospin limit can be found in 
Ref. [32j, however, no analytic formulae for s-wave scat- 
tering lengths were explicitly provided. We, therefore, 
will present it for isospin 1=1 channel. To this end, 



3 Zhi-Hui Guo kindly give the following formulae to me to fitting 
our lattice data. In there, we especially thank him. Without his 
kind help, we can not finish this work smoothly. Nevertheless, 
in this work, we use f n = l32MeV instead of 93 MeV in Guo's 
original formula. 
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-j r _ a 2 \ 32m^- /Tig"! clearly depend on the chiral scale u with chiral logarithm 

+ j^ [s _± mK) ——. ( B3 ) tcrms , 

And the scattering length can be written in a compact 

form, 1=1/ \ q 2 , m K , m £ 2m|- + Sm 2 m 2 

r=i m K f 16m,, 32771^- ^ 9m 4 _ 16m 2 m 2 _ 56m 4 4m 2 _ m 2 

°** - ^l 1+ ^r L5(A1)+ T WA1) + i8«-^) log ^-^ 

'=V,,A fR/n 4 , - , , 9N 256?r 2 m 2 



») > (B4) _^ + 7 7rr ,( s = 4m 2 ,)i^^. (B5) 



(47r/7r) 2^v- T v / -- 

where x^^(^) is the known functions at NLO which 
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